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In this paper, we derive the beta-integral 
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(I could not find this integral in literature). We construct the system of orthogonal 
polynomials related to this integral and also systems related to the Askey integral 
(0.5) and the 5 i? 5 -Dougall formula (0.6). It turns out that all these systems are 
the Wilson polynomials (see [Wil] , [AAR] , [KS] ) outside the domain of positivity 
of the classical weight. 

0.1. Beta-integrals. The integral (0.1) is a representative of a large family of 
so-called beta-integrals. We mention only the beta-integrals that in either case will 
appear in this paper: the de Branges-Wilson integral [dB], [Wil] (see also [AAR], 
3.6) 
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the Second Barnes Lemma (see [AAR], Theorem 2.4.3) 
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the Nassrallah-Rahman integral (see [Ask2], [RS]) 
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and the Askey integral (see [RS]) 

J_ f°° r(l - 2s)r(l + 2s) rfg = £(oi + a 2 + a 3 + a 4 - 3) 

In the last case, the integrand has poles on the integration contour at the points 
s = ±n/2, n = 1,2,.... We can understand this integral in two ways. 

The first way. Fix a such that < a < 1. Consider the integral from — M — a 
to N + a, where M, N are positive integers. We understand this integral as the 
principal value near each pole. For Re a j < 3, the limit of this integral as 
M,N -> +oo exists. 

The second way. Consider the integral from — M — a to N + a over the contour 
that passes the poles from above and consider the limit as M, ./V — > +oo. The both 
variants coincide since the integrand is an even function (the corrections in the 
pairs ±k/2 of poles cancel). 

More complete lists of beta-integrals are contained in [Ask2], [RS], there are 
also numerous g-analogs (see [AW], [Ask2], [RS]), discrete analogs and multivariate 
analogs (see [Gusl], [Gus2]) of beta-integrals. 

Among discrete analogs, we need the Dougall formula 

a + n sin 2na T(ai + a 2 + a 3 + a 4 — 3) 



E 



Jill T ( a J + a + n ) r ( a J -a-n) 2n IliC7<fc<4 r ( a J + a k ~ 1) ' 



Remark. In formulas (0.1)-(0.4), we assume that the parameters aj, bj are real 
and positive. To be definite, consider integral (0.2). Represent it in the form 

1 f°° UUi r K + is)T(a k - is) _ Ui <k <K4 r ( ft fe + «0 



2tt J-oo T(2is)T(—2is) T(ai + a 2 + a 3 + a 4 ) ' 

By the analytic continuation, this equality is valid for Reaj > 0. We can replace 
this restriction by a weaken variant a k + a; ^ 0, —1, —2, . . . (for all the pairs k, I), 
but in this case we must change the contour of integration (see Subsection 0.4) or 
add the sum of residues in the right-hand side of the equation. We can understand 
all the integrals (0.1)-(0.4) and (2.2)-(2.5) in a similar way. 

0.2. Nassrallah Rahman integral. Our derivation of the integral (0.1) is very 
simple, but the integral itself can be obtained by degenerations and hypergeometric 
transformations from the following wonderful Nassrallah-Rahman integral 
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T(a + f)T(f - a)T(g + a)T(g - a) 



ds 



x 



r(l + a - b)T(l + a - c)T(l + a - d)T(X + a)T(X - a)T{p + a)T(jj, - a) 



2a, 1 + a, a + b, a + c, a + d, a + f, a + g, 1 + a — A, 1 + a — /j, 
a, 1 + a — b, 1 + a — c, 1 + a — d, 1 + a — /, 1 + a — g, a + A, a + /i ' 
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+ ..., 

where • • • denotes the sum of the similar summands, where a changes by b, c, d. 

We discuss several integrals (2.2)-(2.5) intermediate between (0.2), (0.1) and 
the Nassrallah-Rahman integral. All these integrals can be obtained from (0.7) 
by degenerations and hypergeometric transforms, however I never have seen the 
final formulas (2.2)-(2.5). Nevertheless they deserve to be written (especially the 
integral representation (2.2) for 3F 2 (1)), furthermore, our calculations are very 
simple conceptually and technically. 

0.3. Finite systems of orthogonal polynomials. We write explicitly the 
system of polynomials p n (s 2 ) orthogonal with respect to the weight 



toi(s) 



2 

IL=i r ( a fc +is) 



T(2is)T(b + is) 



ds. 



Observe that this weight decreases as 



-2(6— ai -02-03)-! 
* 1 



hence only a finite number of the moments 

1 

s 2k wi{s)ds 



f 

Jo 



exists. Thus our system of orthogonal polynomials is finite by the definition. 

We also construct the system of polynomials orthogonal with respect to the 
discrete weight 

where 8 n is the delta-function supported by the point n, and c„ are the summands 
of the Dougall formula (0.6). 

This weight also polynomially decreases. For some aj and a the weight is pos- 
itive, in some cases it is sign indefinite. The latter is not very important for play 
formulas. 

Finally, we construct a system of polynomials (again, this system is finite) associ- 
ated with the weight Ws(s), where Ws(s) is the integrand in the Askey integral (0.5). 
This weight is not positive; moreover, it is not a measure, but we can integrate with 
respect to it. 

It turns out that all the three systems consist of the Wilson polynomials. Fur- 
thermore, the systems related to the weights w 2 (s) and Ws(s) coincide (up to a 
correction factor in the orthogonality relations) . 

In first time, finite systems of orthogonal polynomials (with respect to a con- 
tinuous weight) were discovered by Romanovski [Rom] in 1929 (he constructed 
analogues of the Jacobi polynomials). Askey [Askl] obtained the system of polyno- 
mials related to the Ramanujan beta- integral. Various finite families of orthogonal 
polynomials were considered in a series of papers of Lesky [Lesl]-[Les4]. Some 
applications of such constructions are contained in [Pee] , [BO] . 

0.4. Notation. (a)k ■= a(a + 1) . . . (a + k — 1) is the Pochhammcr symbol. 



q -fp 



ai,.. .a q _ 
61,..., bp' 
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is the generalized hypergeometric function. Next, 

r( Ul )...r(u,) 



Ml, • ■ -U[ 

Vi,...,V m 



T( Vl ) . . .r( Vm )- 

Mellin-Barnes integrals are denned by 



2iri 7_ ; 



ai + s, . . . , a m + s, bi - s, . . . , b n - s 
ci - s, . . . , Cfe - s,di + s, . . . ,di + s 



z s ds. 



We suppose that the integration is given over the contour passing from —zoo to 
+ioo and separating the left series of the poles 



J, 



where j = 0, 1, 2, . . . , oo, 



s = -ai — j, . . . , —a 

from the right series of the poles 

s = b 1 +j,...,s = b n + j. 

We assume that none of poles from the left scries coincides with a pole from the 
right series (this implies the existence of a contour) . Using the standard rules, we 
can represent a Mellin-Barnes integral as a finite sum of hypergeometric functions 
with T-factors. This rules (for m + n ^ k + I, bellow this is always satisfied) are 
contained in the books of Slater [Sla] and Marichev [Mar] , see also [PBM] . 

§1. Preliminaries and preliminary calculations 

1.1. Index hypergeometric transform. Fix a, b > 0. Let / be a function 
on the half- line x ^ 0. The index hypergeometric transform Jb, c is defined by the 
formula 



(1 



•1) J/(.) = ^^~ /(«),*(«. 



is, a- is;a + b;-x)x a+b 1 (1 + x) a °dx. 



\a — b r 



The inverse transform is given by 
(1.2) 

1 



J-^g{x) = 



g{s)2F\{a + is, a — is;a + b; —x) 



TrT(a + b) J 
The transformation J is a unitary operator 

£ 2 (R + ,a; a + b - 1 (l + 2;) a - 6 ) -> L 2 (l 



T(a + is)T(b + is) 



T(2is) 



ds. 



T(a + is)T(b + is) 



T(2is) 



The unitarity condition (the Plancherel formula) in the explicit form is 
(1.3) 

T(a + is)T(b + is) 
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Jo 



fi{x)f 2 {x)x 



a+b-1 



{l + x) a -"dx 



n Jo 



9i(s)92{s) 



T(2is) 



ds 



On properties of this transform, see [FJK], [Kool]-[Koo3], [Yak], [Ner]. The opera- 
tor J is called also the Olevsky transform, the Jacobi transform, the Fourier- Jacobi 
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transform and the generalized Fourier transform, it was discovered by Hermann 
Weyl in 1910. 

1.2. Gauss hypergeometric function. Wc use some simple properties of 
the hypergeometric functions 2 Fi, in particular: the Gauss formula (see [HTF1], 
(2.1.14), [AAR], Theorem 2.2,2) 

, s , s T(c)T(c- a-b) 

(1.4) 2 F 1 (a,b;c;l) = r c _ o)r(c _ b i Re(c-a-b)>0. 
and the Bolza formula (see [HTF1] (2.1.22-23)) 

(1.5) 2 F 1 (a,b;c;-x) = (1 + x)- a 2 F 1 (a, c - b; c; -?—) = 

x + 1 

(1.6) = (l + .T)- b 2 F 1 ( C - a ,6;c;-^) = (1 + x) c - a -\F 1 (c - a, c - b; c; -x) . 

Recall that the function 2 Fi(a,b;c;x) admits the analytical continuation to the 
half-line x < 0. The asymptotics of the function 2 F\(a + is, a — is; a + b; —x) as 
x — > +00 has the form 

2 F\(a + is, a — is; a + b; —x) = ip(x)x~ a , 

where ip is a bounded function for s ^ and ip(x)/ \n(x) is a bounded function for 
s = (for a more explicit expression, see [HTF1],(2.3.2.9)). 

The asymptotics of the function 2 i*i (a + is, a — is; a + b; —x) as s — > +00 has the 
form 

2 Fi(a + is, a - is; a + b; -x) = \{s) S - a - b+1/2 . 

where A(s) is a bounded function (see the Watson formula [HTF1], (2.3.2.17)). 

Recall three representations of the Gauss hypergeometric function as Barnes 
integrals see [HTF1], (2.1.15), [Mar], §2.10, formulas 11(1), 19(1), 22(1), [AAR], 
Theorems 2.4.1, 2.4.2, [PBM], 8.4.49(13), (20), (22); below x > 



2tti r(a)r(6) 7_ loc r(c-s) 



(1.8) 

1 r(c) 



2tti r(a)r(6)r( c 



^— — — - / r( s )r( s + c - a -6)r(a- s )r(6- s ) a; - s d s = 

- a)r(c- 6) 7_ loo 

= 2 Fi(a,b; c; 1 - x); 



(19) 1 F(c) r r Wa + c-a-6) 

(i - 9J 2« i(cj l 40 or( s + c -a)r( s + c -6) a: ds ~ 



_ J(l-x) c - 1 2 Fi(a,6;c;l-x); x > 1, 
0; x < 1 
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1.3. T-function. We shall use the following integral representation of the beta- 
function 

r°° x^dx 

B(a,p-a), 



I 

Jo 



o (l+^) p 

The substitution x/(x + 1) = z transforms this integral to the usual definition of 
the beta-function. 

Asymptotics of the function T(z) is given by the Stirling formula 

z — > oo arg z < ir — e. 



T(z) = V2^z z - 1/2 e- z (l + 0(1/*)); 
In particular (see [HTF1], (1.18.3), (1.18.6)), 



r ( Z + fl ) a-b 

~ z 

T{z + b) 

|r(a + *s)| 



2n\s 



oo, 



o-l/2 e -kK/2 



arg z < 7r - e, 



1.4. Mellin transform. Consider a function / on the half-line x ^ 0. Its 
Mcllin transform is defined by the formula 



*(*) = / 
Jo 



x s f(x) 



dx 



The inversion formula is 



The Mellin transform takes the convolution 

r°° dt 

J o h(x/t)f 2 (t) T 

to the product of the Mellin transforms /i(s)/2(s). For precise statements, see for 
instance [Mar]. 

Let a > 0. Obviously, the function f(ax) is mapped to a~ s f*(s). Also, the 
function x a f(x) corresponds to f*(s + a). 
1.5. We need some integrals of the form 



(1.10) 



o {x + z)p 

1 

2^1 



2 F 1 (p,q;r;-x) dx = 





r 


pioo 

r 






J —too 



s + a, p — s — a,p + s,q + s, —s 
r + s 



z s ds 



The Mellin-Barnes integral in the right-hand side can be expressed as a linear 
combination of two 3 F 2 . 

As a quite standard (see [Mar]) pattern, we prove the identity (1.10); this state- 
ment also is quite standard (see [PBM], formula (2.21.1.16)). For this, represent 
x a (x + z)~ p as the inverse Mellin transform 



x a {x + z)- p = 



?a-p 



2nr(p) 



{ 



T(s + a)T(p -s- a)z s x s ds 
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(for a verification, it is sufficient to evaluate the direct Mellin transform). Applying 
the Barnes integral (1.7) and the convolution formula for the Mellin transform, we 
obtain (1.10). 

1.6. Corollaries from formula (1.10). 

Lemma 1.1. a) The operator J takes the function 



(l + z)-°- 



to 



T(c + is)T(c-is) 



T(c + a)T(c + b) 
b) The operator J takes the function 

(l + x) b - a T(c + is)T(c-is) 



to 



2F 1 



c + is, c — is 
c + a 



(x + z) c + b "' T(c+a)T(c + b) 

c) The operator J takes the function 

x ~u-a t() r(-U + b) r(u + l s)T(u- iS ) 

T(a + u) T(b + is)T(b-is)' 



Proof. We must substitute these / to formula (1.1). We can simply refer to the 
rows (2.21.1.4), (2.21.1.5), (2.21.1.15) of Prudnikov, Brychkov, Marichev, vol.3, or 
to the rows (7.512) of Ryzhik-Gradshtcyn. 

But it is better to explain what happened and why for some values of the pa- 
rameters the integral (1-10) admits an explicit evaluation and for some values this 
is impossible. 

1. For a = r in the right-hand side of (1.10), two T-factors in the integrand 
cancel. The remains is a Barnes integral of the type (1.8); this proves the statement 
b). 

2. Substituting z = 1 to the statement a), we obtain the claim b). 

3. Substitute z = 1 to (1.10). For r = p + q + p, the right-hand side can be 
evaluated by the Second Barnes Lemma (0.3). This gives c). H 

1.7. Integrals with products of hypergeometric functions. A direct 
application of the convolution formula for the Mellin transform gives 



(1.11) 



/ x a ~\Fx 


p,q 

; —lux 


2F1 


u,v 

■ — cox 


dx 


h 


r 




w 





27ri W 



r, w 


pioo 




r 


u,v,p,q_ 


J —zoo 



a + s,u + s,v + s,p — a — s,q — a 
r — a — s,w + s 



ds. 



Substitute a = w = r to this identity. Then four T-factors of the integrand in the 
right-hand side cancel. Applying the Barnes formula (1.8), we obtain 



(1.12) 



r x t - i 2 f 1 


p,q 

: —ojx 


2F1 


u.v 

■ —LUX 


dx 


lo 


r 




r 





= uo- r+u Cj- u T 



r,r,p — r + u,q — r + u,p — r + v,q — r + v 
u, v,p, q,p + q + u + v — 2r 



x 2 F! 



p — r + u,q — r + u _ 
p + q + u + v — 2r '' 
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It is formula (2.21.9.7) from Tables [PBM] (where a prime in the second index of 
hypergeometric function in the right-hand side is lost). 
We also need the integral 



(1.13) f/^fl-^-Vi 
Jo 









V 







dz = 



= — r 

2ni 



jj, + s, n + v — a — P + s,(f> + s,ip + s, — s, £ — — ip + s 
v + ii — a + s, /i + v — (3 + s 



ds. 



To deduce this formula, we represent (1 — z) v ~ 1 2F 1 {a, (3; v; 1 — z) as the integral 
(1.9), and represent 2-Fi(</>, ip: 1 — z ) as the integral (1.8). It remains to apply the 
convolution formula. 

1.9. Corollaries from the formula (1.11). 

Lemma 1.2. a) The transform J takes 



p + b,q + b 
a + b ' 



to 



y b -*T{a + b) 



(l + xf- a 
T(p + is)T(p - is)T(q + is)T(q - is) 



T(p + q )r(p + b)r( q + b) 

b) The transform J takes 2 F\ 
T{a + b) 



T(a + is)T(a — is) 
p + b,q + b 



(1 + x) 



b—a 



Fi 



to 



p + is,p — is i 

p + q 



a + b ' 

T(p + is)T(p — is)T(q + is)T(q — is) 



T(p + q)T(p + b)T(q + b) T(a + is)T{a-is) 

c) The transform J takes 2-F1 a + c i a + ^ x ^ Q 

[a + b + c + a 

r(a + b + c + d) ■ T(c + is)T(c - is)T(d + is)T(d - is) 
r(a + c)r(a + d)T(b + c)T(b + d)T(c + d) ' 

Proof. The claim a) is an immediate corollary of formula (1.11). The substitution 
y = 1 to a) gives the statement b) (of course, it is more pleasant to substitute ui = Co 
directly to (1.11). 

Let us deduce c). Substituting / = 2-Fi[---] to (1.1) and applying the Bolza 
transformation, we obtain 



r( 



1 r°° 

4tm / ^^Fi 
a + b) J 



b + d,b + c 
a + b + c + d 



\-x 



a + is, a — is 



a + b 



;-x 



dx 



Furthermore, in the integrand in the right-hand side of (1.12), two T-factors cancel 
and we obtain 



J-l 

2ni 



a + b + c + d 
a + is, a — is,b + c,b + d 



/ 1 r 

-ioo 



a + is + t,a — is + t,d — a — t,c — a — t,—t 
c + d-t 



dt 



BETA-INTEGRALS AND FINITE ORTHOGONAL SYSTEMS OF WILSON POLYNOMIALS9 



It remains to apply the Second Barnes Lemma (0.3). 

§2. Degenerate cases of the Nassrallah Rahman integral 

In this section, we apply the Planchcrel formula (1.3) for the index hypergeo- 
mctric transform to various pairs of functions from Lemmas 1.1-1.2. 

To be brief, we assume that the parameters a, b, c, d, e, /, u, v are real positive. 

2.1. De Branges Wilson integral. Apply the Plancherel formula (1.3) to 
the pair of functions (1 + x)~ a ~ c and (l + x)~ a ~ d , see Lemma 1.1. In the left-hand 
side we have the beta-integral 



(2.1 .) 



roc x a+b- 
J (l + x)** 



roo ^+b-l dx 

i+b+c+d 



In the right-hand side, we obtain the de Branges- Wilson integral (up to T- factors). 
Evaluating (2.1), we obtain (0.2). 

This proof of the de Branges- Wilson integral is known, see the work of Koorn- 
winder [Koo2]. 

2.2. The integral (0.1). It is sufficient to apply the Plancherel formula to the 
pair of functions x~ u ~ a and x~ v ~ a , see Lemma l.l.c. 

2.3. Apply the Plancherel formula (1.3) to the pair of functions 



(l + x)- a - e and 2 F X 
In the left-hand side we have 



a + b + c + d' 



/ x a+b - 1 {l + x)- b - e 2 F 1 
Jo 



a + c, a + d 
a + b + c + d'' 



dx. 



We transform this by the Bolza formula and obtain 



a + c,b + c 
a+b+c+d 



dz. 



This is a standard integral representation for 3 F 2 (see [AAR], (2.2.4)), and finally 
the left-hand side is 



T(a + b)T(c + e) 
T(a + 6 + c + e) 



3-^2 



a + c,b + c,a + b 
a + b + c + d, a + b + c + e' 



Equating the left-hand and right-hand sides of the Plancherel formula, we obtain 



( 2 - 2 \ I 
n Jo 



T(a + is)T(b + is)T(c + is)T(d + is)T(e + is) 



ds 



T(2is) 

T(a + b)T(a + c)T(a + d)T(a + e)T(b + c)T(b + d)T(b + e)T(c + d)T(c + e) 
T(a + b + c + d)T(a + b + c + e) 

a + c,b + c, a + b 



x 3 F 2 



a + b + c + d, a + b + c + e' 
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Among the T-factors in the numerator, the factor T(e + d) is absent. 

Remark. Obviously, the left-hand side of the identity is symmetric in a, b, c, d, e. 
Symmetry of the right-hand side is equivalent to the Kummer formula, see [AAR] , 
Corollary 3.3.5. 

2.4. Now apply the Plancherel formula (1.3) to the pair of functions 



2*1 



a + c, a + d 



, ; —x 



and 



a + b + c + d ' 
see Lemma 1.2.c. In the left-hand side, we have 



*i 



a + e, a + f 
a + b + e + f 



f 

Jo 



x 



a+b—l 7-1 
2-f 1 



-f 

JO 



a + c, a + d 
a + b + c + d 

c a+b-i 2Fi 



;-x 



2*1 



a + e,a + f 
a+b+e+f 



; -x 



(l + x) a -"dx = 



b + d,b + c 
a + b + c + d'~ 



2*\ 



a + e,a + f 



r 



dx. 



Applying formula (1.11), we obtain 



(2.3) 



1 



T(a + is)T(b + is)T(c + is)T(d + is)T(e + is)T(f + is) 



ds 



T(2is) 

^-T(a + c)T(a + d)T(c + d)T(b + e)T(b + f)T(e + f) x 

a + b + s,a + e + s,a + f + s,d — a — s,c— a — s,—s 
c + d— s, a + b + e + f + s 



ds. 



The expression in the right-hand side can be represented as a linear combination 
of three functions 4*3(1). 

2.5. Finally, we apply the Plancherel formula to the pair of functions 



2*1 



p + b,q + b 



(1 + x) 



b—a 



and 2*1 



u + b, v + b 



; -x 



(l+x) 



b—a 



(their index transforms are evaluated in Lemma 1.2b). 
We must evaluate the integral 



[ 

Jo 



- a+b - 1 (l+x) b - a 2 F 1 



p + b,q + b 
a + b ' 



2*1 



u + b, v + b 



a + b 



; -x 



dx. 



Using the Bolza formula, we obtain 



[ 1 y a+b - 1 (l-yr +p 2 F 1 
Jo 

-f 

Jo 



p + b,a- q 
a + b 



\y 



(l-z) 



a+b— 1 u+p 



p + b,a- q 
a + b 



2*\ 



;1 



u + 6, a — v 



;y 



dy = 



2^1 



u + b, a — v 
a + b 



;i 



dz. 



It remains to apply (1.13). As a result, we get 



(2.4 .) 



r(6 + is)T(p + is)T(q + is)T(u + is)T(v + is) 



2iri 



T{2is)T(a + is) 

u + v,p + q,p + b,q + b 

a — v,u — v 



ds 



f 



u + p + s,u + q + s,b + u + s,a — v + s,v — u — s,—s 
u + a + s,u + b + p + q + s 
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Remark. The right-hand side can be written in the form 

r [u + v,p + q,p + b,q + b~\ x 

a — v,u — v 

u + p,u + q,u + b,a — v 



v 



u,u + p,u + q,u + b,a 
u + a,u + b + p + q 



4F3 



+ r 



u — v,p + v,q + v,b + v,a — u 
v + a, v + b + p + q 



4F3 



1 + u — v,u + a,u + b + p + q' 

p + v,q + v,b + v,a — u 
1 — v + u,v + a, v + b + p + q' 



+ 



Rahman [Rah] gives the for integral (2.4) a 7_F 6 -expression (see also [GR], (6.3.11)). 
However, it can be easily reduced to the form 4F3 using the nonterminating rF e - 
Whipple transform, discussed in Bailey's book [Bai], 4.4, 6.3, 7.5. 

Remark. The left-hand side of the identity (2.4) is symmetric in b,p,q,u,v. In 
the right-hand side, the symmetry in u, v and in b,p, q is obvious. The transposition 
of b and u gives a 4-term .^-identity. Substitute 

a = v — to, 

where m is a nonnegative integer, to this identity. Then two summands vanish (due 
r(— m) in the denominators), and we obtain 4F3-Whipple transform (see [AAR], 
Theorem 3.3.3; it is mentioned below in Subsection 3.1). 

2.6. Nassrallah Rahman integral (0.4). If 

a = b + u + v+ p + q, 

then two T-factors in the right-hand side of (2.4) cancel. Applying the Second 
Barnes Lemma (0.3), we obtain (0.4). 

2.7. Multiplying the both sides of (2.4) by T -2 (&) and passing to the limit as 
b — ► +00, we obtain 



(2 



1 r 

n Jo 



T{p + is)T(q + is)T(u + is)T{v + is) 



T(2is)T(a + is) 



ds = 



2iri 



u + v,p + q 

a — v.u — v 



u + p + s,u + q + s,a — v + s,v 
u + a + s 



ds 



(in the right-hand side we have a linear combination of two functions 3^2(1)). 
§3. Finite systems of orthogonal polynomials. 

3.1. Wilson polynomials. In the famous work [Wil], Wilson constructed the 
polynomials p n {s 2 ) orthogonal with respect to the weight 



w(s) = — 

7T 



they are defined by the formula 
(3.1 .) p n (a,b,c,d;s 2 ) = 

= {a + b) n (a + c) n (a + d) n 4-F3 



r(o + is)T(b + is)T(c + is)T(d + is) 



T(2is) 



-n,n + a + b + c + d— 1, o + is, a — is _ 
a + 6, a + c, a + d 
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The orthogonality relations have the form 

poo 

(3.2) / p n (a,b,c,d;s 2 )p m (a,b,c,d;s 2 )w(s)ds = 
Jo 

_ nW(a + b + n)T(a + c + n)T(a + d + n)T(b + c + n)T(b + d + n)T(c + d + n) 
~ T(a + b + c + d + n)(a + b + c + d + 2n - 1) m '' 

Evidently, the polynomials p n (a, b, c, d; s 2 ) are symmetric with respect to b, c, d. The 
symmetry in four indices a, b, c, d is equivalent to 4i 7 3-Whipple transform mentioned 
above in Subsection 2.5. 

3.2. Proof of orthogonality. Several proofs of orthogonality of the Wilson 
polynomials are known. Our purpose is to give a model proof that works in 3 cases 
of " exotic orthogonality" discussed below. 

Lemma 3.1. Let V , W be N -dimensional linear spaces with bases €j, fj respec- 
tively. Define the sesquilinear form B(-, •) on V x W by 

(3 ' 3) B{ekJl) - f{^TVTkTTy 

where [i, v are fixed. Then the system of vectors 

(-n) 3 (n + fi + v -1)3 ^ _ ^ {-n)j{ji + n + v- l) 3 

is biorthogonal, i.e., B{Rk 1 Ti) = for k^l. Moreover, 

n!ro)ro + n) 



(3.4) B(T n ,R n ) r(l/ + n))r(n + Ai + 1/ _ 1)(/i + 1/ + 2 „_ 1) - 

Proof. This lemma imitates the orthogonality relations for the Jacobi polynomi- 
als. Let V — W be the space of polynomials on the segment [0, 1] with the scalar 
product 



1 /-i 

<p,q>= T(pj J Q v( x )^ x ) xV i l -xY <&•■ 



Let e n — f n be the function x n . Then < fi > coincides with (3.3). Next, 
Rn = T n are the Jacobi polynomials P tl ~ l v ~ 1 (2x — 1) in the standard notation 
(see [HTF2], 10.8). □ 
Now let V — W be the space L 2 with respect to the Wilson weight w(s). Assume 

,„ c x ,s (a + is) k {a-is)k , (b + is) k (b - is) k 

(3.5) e k (s) = — — ; — , fk(s) - 



T(a + c + k)T{a + d + ky J y ' T(b + c+k)T(b + d+k)' 
Then by (0.2), 

/ e k (s)f m (s)w(s)ds = 
Jo 

1 



X 



7iT(a + c + k)T(a + d + k)T(b + c + m)T{b + d + m) 
r(a + k + is)T(b + m + is)T(c + is)T(d + is) " 







ds = 



T(2is) 

_ T(c + d)T(a + b + k + m) 
~ T(a + b + c + d + k + m)' 
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We obtain the relations of the type (3.3) for scalar products. In notation of Lemma 
3.1, we have 

R n = ot n Pn(a, b, c, d; s 2 ); T n = (3 n p n (b, a, c, d; s 2 ), 

where a n ,(3 n are normalizing factors (we omit them to be brief). It remains to 
recall that the Wilson polynomials are symmetric with respect to a, 6, c, d. Hence 
R n and T n coincide up to a factor (precisely this place of the proof is surprising). 
Formula (3.2) follows from (3.4). 

3.3. Orthogonality relations associated with the integral (0.1). Consider 
the weight 



(3.6) w 1 (s) = - 

7T 



on the half-line s ^ 0. 
Assume 



T(p + is)T(u + is)T(v + is) 



T(2is)T(q + is) 



2 



Then 



Jo 



( \ r(-7i -k + q) 

e "( s ) = r( P + u + fc) (u + ts ' k{ - u ~ ls ' k] 

fn{8)= T( P + v + k) (» + ")*(»-»)*■ 



,,,,, , v, T(u + v + k + l)T(q - u - v -p - k - I) 
e fc (s)/ i (s)«;i(s)ds = ^ _ ^ = 

(-l) k+l ir T(u + v + k + l) 



T(q — p) sinir(q — u — v — p) r(l — q + u + v + p + k + I)' 



Observe that the scalar products have the form (3.3). Thus we can repeat literally 
all the considerations of Subsection 3.2. 

There exists a way that is even more simple. 

Lemma 3.2. Let £ be a linear functional on the space of even polynomials 

c + cis 2 H h c N s 2N 

and for 

hk = (a + is) k (a- is) k 

we have 

T(a + b + k)T(a + c+k)T(a + d + k) 
(3J) l{hk) = C T(a + b + c + d + k) • 

Then for k + I < N the Wilson polynomials Pk,Pi satisfy 

C 

£(p k (a, b, c, d; -)pi(a, b, c, d; ■)) = — — — — — — — • a k (a, b, c, d)S k ,i, 

l(o + c)l (o + d)l (c + d) 
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where a n (a,b,c,d) is the expression situated in the right-hand side of (3.2). 

In other words, the Wilson polynomials form an orthogonal system with respect 
to the scalar product 

(f,g)=e(m 

in the space of polynomials. 
Proof. Assume 

roc 

= / g(s)w(s)ds, 

Jo 

where w(s) is the Wilson weight. Then (3.7) is satisfied and our statement is a 
rephrasing of the orthogonality relations for the Wilson polynomials. □ 

Now suppose 

roc 

= / g(s)wi(s)ds 
Jo 

and 

hk(s) = (p + is) k (p-is) k . 

We have 

, . _ r(p + v + k)T(p + u + k)T(u + v)T(q -u-v-p-k) _ 
( k >~ T{q-p-k)T(q-u)T(p-v) " 

T(u + v) sin Tr(p - q) T(p + v + k)T(p + u + k)T(l -q + p + k) 



T(u — q)T(q — v) sin ir(u + v + p — q) T(l — q + u + v+ p + k) 

Therefore, the system of the Wilson polynomials 

p n (p,u, v, 1 - q; s 2 ); An<q-p-u-v-l 

is orthogonal with respect to the weight (3.6). 

3.4. Orthogonality relations associated with Dougall formula. Consider 
the weight 

oo ot + t 
W2(t) = —i S(t - n), 

where S(s — n) is the delta-function supported by the point n. For uniformity, 
substitute t = is — a and transform this expression to the form 
(3.8) 

4 . / , \ oo , 4 

It _.[ 

j=l " -oo k j=l 



4 • ( 4- \ °° ( 4 1 

^ _ = j-j smn^ + a) _ ( . s) ^ J -Q r(1 _ a . _ ls)r{1 _ a . + ls) . 5{ls L 

i' = 1 — on ^ i = 1 ' 



To apply Lemma 3.2, suppose 

t(9) = J g(s)w 2 (s)ds 
h k {s) = (1 - oi + is)fe(l - ai - is) k . 
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Then 

f 39 ) e(h .)- sin(2 TO ) x 

( ' 1 k) ~ 2^r( fl2 + a 3 )i> 2 + a 4 )I> 3 + 04) 

T(a 1 + a 2 + as + 0,4 — k — 3) _ 
r(fl! + a 2 - k - l)r(cn + a 3 - k - l)T(ai + a 4 - fc - 1) ~~ 

sin(27ra) sin 7r(ai + a 2 ) sin n(ai + a 3 ) sin 7r(ai + a 4 ) 
27rr(a 2 + a 3 )r(a 2 + a 4 )r(a 3 + a 4 ) sin7r(ai + a 2 + a 3 + a 4 ) 

r(-ai - a 2 + k + 2)r(-ai - a 3 + k + 2)r(-ai - a 4 + k + 2) 
r(4 — ai — a 2 — a 3 — a 4 + fc) 

Thus we obtain that the system of the Wilson polynomials 

(3.10) - a 4 , 1 - a 2 , 1 - a 3 , 1 - a 4 ; s 2 ); 4n< a\ + a 2 + a 3 + a 4 - 3 

is orthogonal with respect to the Dougall weight (3.8). 

3.5. Orthogonality relations associated with the Askey integral (0.5). 

Now consider the Askey weight, i.e., the weight that is the integrand in (0.5) 

r(i-2s)r(i + 2s) 



Il-=i T( aj + s)T( aj - s) 



In notation of Lemma 3.2, 

/oo 
g(s)w 3 (s)ds, 
-OO 

h k (s) = (l-ai + s)k(l - 01 - s) k . 



Then 

£(h k ) = 



1 



r(a 2 + a 3 )T(a 2 + a 4 )r(a 3 + a 4 ) 

r(ai + a 2 + a 3 + a 4 — k — 3) 



r(ai + a 2 - k - l)r(ai + a 3 - fc - l)r(a 4 + a 4 - fc - 1) 



and we obtain the expression coinciding with (3.9) up to a constant factor. Hence, 
we obtain the same system (3.10) of orthogonal polynomials. 

§4. Examples of index integrals 

Index integrals are by themselves fairly known subject, see, for instance [Sla]. 
In this section, we present several amusing index integrals extending the integrals 
(0.1), (0.2), (2.2)-(2.5). 

4.1. The index transform of the function (1 + x) b ~ a (x + y + l)~ c ~ b was evaluated 
in Lemma 1.1. Using the inversion formula, we get 



T(a + is)T(b + is)T(c + is) 



T(2is) 



iFi 



c + is,c — IS _ 
a + c 



iFi 



a + b 



TrT(a + b)T(a + c)T(b + c) 
(l + x + y) c+b ' 
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4.2. Writing out the Plancherel formula for the functions (l+x) b a (x+y+l) c b 
and (1 + x)~ a ~ d , we obtain (see de Branges [dB], Theorem 13) 



T(a + is)T(b + is)T(c + is)T(d + is) 



T(2is) 



c + is.c — is 



a + c 



-y 



ds 



TrT(a + b)T(a + c)T(a + d)T(b + c)T(b + d)T(c + d) 
T{a + b + c + d) ; 



Fi 



b + c, c + d 
a + b + c + d 



\-y 



This formula coincides (up to a permutation of letters) with the inversion for- 

see Lemma 1.2.c (this implies the Second Barnes 



mula for 2F1 



a + b + c + d 



;-x 



Lemma (0.3); evidently, it is not the most simple its proof). 

4.3. Applying the Plancherel formula to the pair of functions (1 + x) b ~ a (l + x + 
y)- b - c and (1 + x) b - a (l + x + z)- b - d , wc obtain 



n Jo 



T(a + is)T(b + is)T(c + is)T(d + is) 



T(2is) 
x 2-F1 



c — is, c + is 
a + c 



;-y 



2F1 



d — is, d + is 



T(c + a)T{c + b)T(d + a)T(d + b) 



a + d 
x a+b - 1 (l + x) b - a dx 



ds = 



l (x + y + l) b+c (x + z + l) b + d ' 

In the right-hand side we have one of integral representations of the Appel function 
F\. Substituting z = y, we obtain 



t Jo 



T(a + is)T(b + is)T(c + is)T(d + is) 



T(2is) 

X 2^1 



c — is,c + is 
a + c 



\-y 



iF\ 



d — is, d + is 
a + d 



ds = 



7rr(a + b)r(a + c)T(a + d)T(b + c)T(b + d)T(c + d) 
_______ 



-2F1 



2b + c + d,c + d 
a + b + c + d '■ 



4.4. Applying the Plancherel formula to the pair of functions (l+x) b a (l + a; + 



y)- e - b and 2 Fx 



a + c,a + d 
a + b + c + d' 



we get 



T Jo 



T(a + is)T(b + is)T(c + is)T(d + is)T(e + is) 



T(2is) 



2F1 



a + is, a — is _ 
a + c 



ds = 



2ni 



7 (1 + y) a ~ e T(b + c)T(b + d)T(c + d)x 



- 1 r 



J 



a + b + s, a + c + s,a + d + s,e — a — s, —s 
a+b+c+d+s 



(i + y) s ds 



(it is a sum of two functions 3-F2). 

4.5. It is easy to extend this list. Even our Lemmas 1.1-1.2 were not completely 
utilized. 
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